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By F. N. H. ROBINSON 

(Manuscript received January 18, 1967) 

We consider, from a number of different viewpoints, the tensor coefficients 
which describe second harmonic generation, optical rectification, and the 
Pockets or linear electro-optic effect in acentric crystals. Stationary per- 
turbation theory is used to calculate the low-frequency limit of the intrinsic 
electronic nonlinearity neglecting all effects due to local fields or lattice 
polarization. Solid methane is used as an example and the result used to 
estimate the coefficient in hexamethylenc tetramine. The calcidated result 
is within a factor of 2 of the experimental figure. The method is susceptible 
to further refinement and, since it requires only a knowledge of ground 
state wave functions, and is essentially very simple, it appears to offer a 
useful approach to the calculation of the coefficients. 

The classical anharmonic oscillator model is briefly covered and the model 
is related to a quantal treatment. We find that the anharmonic potential 
used in the model is directly related to the actual crystalline potential. It 
can also be related to the charge distribution in the electronic ground state. 

Local field corrections and the effects of lattice polarization are presented. 
These alter the nonlinear properties in a simple and obvious way, but one 
which has been misunderstood in some of the literature. 

Our results form a theoretical background to Miller's empirical ride 
relating the nonlinear coefficients to the linear susceptibilities. An extensive 
table of Miller-reduced tensor coefficients collated from the published litera- 
ture is presented. 

Finally, we draw together some of the threads of the previous sections. 
An appendix deals with the vexing question of definitions. 

I. INTRODUCTION 

Second harmonic generation, optical rectification, and the linear 
electro-optic effect are particular aspects of a process in which two 
fields, E p , e and EJc' yl , generate a polarization 

P°e ial = d^E^ l E\e iyt . (1) 

913 
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Our concern is with the tensor coefficients d?f? which (Nye 1 ) necessarily 
vanish in centric (centrosymmetric) crystals and which, in acentric 
crystals, are subject to symmetry restrictions, (Kleinman 2 ) which often 
leave only one or two independent components. 

Experimentally, the values of the allowed components of d in different 
materials and at different frequencies range from about 2.10 -1 esu 
(cm/stat volt) to about 6.10" 5 esu. This range may be contrasted with 
the linear optical susceptibility x which is between 0.1 and 0.3 for the 
vast majority of materials and only quite exceptionally exceeds unity. 
There is, however, a connection between the tensor d and x which is 
expressed by an important empirical rule due to Miller. 3 If we write 
d-fr as 

tfff - x:*f*iA«* , ( 2 ) 

where x," is the ii component of x at a frequency a, and if we have chosen 
a principal axis system for x, then the allowed components of A,, A 
for all effects and all materials are similar in magnitude. We shall see 
in a later section that for very many materials in both the visible and 
10 /u region of the spectrum (Patel 4 ), A,,* is near 3 X 10" 6 esu. No 
materials with A above 20 X 10~ 6 esu have yet been found and very 
few are known to have a value below 0.2 X 10~ fl esu. In the case of 
NH 4 H 2 P0 4 where the best measurements of s.h.g., optical rectification 
and the electro-optic effect are available (Francois, Ward, 6 Carpenter 7 ) 
the value of A, 23 from all three effects is 3 X 10~ 6 esu within the experi- 
mental error of 15 percent. The fact that s.h.g., a purely optical effect, 
leads to the same value of A as rectification and the electro-optic effect 
indicates quite clearly that the basic mechanism of the nonlinearities 
is common to all three effects and must therefore reside in the electronic 
motion of the system. In the next section, we shall concentrate on this 
aspect of the problem and neglect the effects of local fields and lattice 
polarization. 

A number of authors (see Section IV for references) have given quantal 
treatments of the optical nonlinearities whose end result is an expression 
for the coefficients d'fZ in terms of sums of rather inaccessible matrix 
elements. Useful as these expressions are, in establishing some of the 
general properties of the coefficients, they are not a practical step on the 
road to calculating the coefficients from other empirical quantities. At 
the other extreme, the classical anharmonic oscillator model has been 
used to demonstrate some of the qualitative features of nonlinear be- 
havior (see Section III). This treatment, though simple and appealing, 
suffers from the defect that the relation between the parameters of 
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the model and those of the real system is obscure. In Section IV, we 
shall remedy this defect and show that the two approaches are closely 
related. 

First, however, we give an approximate method of calculating the 
low-frequency limit of the coefficients from stationary perturbation 
theory in a form in which it has been successfully applied to the linear 
properties of n electron atoms (see, e.g., Dalgarno 8 ). 

II. MAGNITUDE OF THE COEFFICIENTS 

At low frequencies, i.e., well below any electronic resonance we can 
use stationary perturbation theory to calculate, to arbitrary order in 
the applied field E the energy W of the ground state. The polarization 
is then given by 

We shall assume that we are dealing with a crystal containing N identical 
atoms or molecules in unit volume whose individual ground state energies 
are w, so that W = Nw. 

If H is the Hamiltonian of an unperturbed molecule, its Hamiltonian 
in the field E is 

H = //„ + h = H n - eER, (4) 

where 

cB = c £ >;"' (5) 

m = l 

is the dipole moment operator of the molecule, and the sum extends over 
all n valence electrons. We can neglect the core electrons because of 
their high binding energies. If we expand w in increasing order in E as 

w «= w -J- Wi + u>2 + w 3 , etc., (6) 

the term w, gives the permanent dipole moment of the molecule, w 2 
gives the linear susceptibility and w 3 gives a polarization quadratic in 
E which leads to the desired nonlinear coefficients. The electric field 
will perturb the state function \f/ and we shall write the perturbed func- 
tion as either 

i> =■ *o + & + * 2 + • ■ • or | > = |0) + |1) + |2) + • • • . (7) 

Knowledge of ^ or | ) to first order in E is sufficient to determine u>, ,w 2 , 
and w% for 
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Wl = (0| h |0> 

w 2 = (0| h |1> (8) 

i*-<l|A|l>. 

Moreover, the correct value of ^ or | 1) is determined by the require- 
ment that it minimize w. Thus, we can obtain | 1) by a variational 
procedure and the only element of choice left to us is that of the trial 
wave function. 

Minimizing w is equivalent (see Dalgarno and Lewis 9 ) to the simpler 
problem of minimizing 

<1|#„ |1> + 2<0| h |1>, (9) 

where the notation fi or h means H - (0\ H„ | 0) or h - (0 | h | 0). 
As a trial function, we take 

|1> = \h |0> (10) 

so that (9) becomes 

X 2 (0| h.HJi |0> + 2X<0| fr |0). (11) 

The minimization with respect to X gives 

_ (o| h 2 |o) m _ (o| k 2 |o) (12) 

<0| hft B h |0) (01 hHh |0> 

The unperturbed Hamiltonian of the system is of the form 



h 

•2) 

and so, in the denominator of (12), 



H = -3T- S V; + r (13) 

2 m „T7f 



£jk-*£. + -aff-EV. • (14) 

m „ 

Thus, 

<0| /*iU |0> = — f *. S fft-^' Z v.*. dr, (15) 

where dr is an element of configuration space and we have used H | 0) = 
0. Equation (15) can be written as 

(oi kAji. |o) = -ir-B' I S v„tf E *•«•• dr, (ie) 
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and integrated by parts, to give 

<0| hltji |0> - +~ B-ElZfti dr, (17) 

where the discarded first integration part vanishes at the limits, if 
these are infinite, or if they are the boundary of a cell in periodic lattice, 
provided only that E does not vary appreciably within a cell (dipole 
approximation) . 

If we are dealing with isolated atoms j \f/ 2 u dr = 1 and we have 

(0| h&Ji |0> = +n ~- E-E, (18) 

a somewhat unfamiliar form of the sum rule. If, on the other hand, 
we are dealing with overlapping molecules in a periodic lattice, the 
variational problem is to minimize the contribution to w from a single 
cell of the lattice. Thus, in (11) and all succeeding equations, the inte- 
grals implied by the expectation values are to be taken only over the 
interior of a cell. This will also apply to all integrals involved in evalua- 
ting w 2 = (1 | h | 0) and w> 3 = (1 | h | 1). In this case (18) remains un- 
changed. This can be shown to be a general consequence of time reversal 
invariance and the commutation rule 

(p, g) = ih. (19) 

We now have 

|D - u |0) = +^& h |0) (20) 



or 



m . + *w$*rn w .* m . (21) 

From this we obtain the second-order energy 

W2 na a E 2 ' (22) 

where a = h 2 /me = 0.53 A. If we let E = E x ,0,0, and R = X, Y, Z 
this gives 

■*- -^-(XJEI (23) 
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and the atomic polarizability is 



a = —{l 2 )\ (24) 

na„ 



For the H atom, this gives a = 4al instead of the correct value 4.5a^, 
while for the helium atom, taking an effective nuclear charge Z = 
27/16 gives 1.8 X 10~ 25 ccs. The experimental value is 2.1 X 10~ 25 ccs. 
In general, (24) is a lower limit to a, if we evaluate (X 2 ) correctly as the 
expectation value of the mean square moment of all the electrons. If 
the electrons are uncorrelated 

(I 2 ) = n(x 2 ), (25) 

where (£ 2 ) refers to one electron. We used this procedure in helium since 
the two electrons are in orthogonal spin states and are automatically 
uncorrelated. In more complicated atoms correlation exists and almost 
always results in 

(X 2 ) < n(x) (20) 

since electrons repel each other. Thus, while (24) is a lower limit we 
cannot say anything about the sign of the error in 

« = ^-{x 2 ) 2 . (27) 

do 

We note, in passing, that, in a solid with overlapping molecules, a the 
polarizability is large. This leads to an element of instability in the situa- 
tion for as a increases the screening of the coulomb potential becomes 
more effective and the electrons less localized leading to a further in- 
crease in a and eventually metallic behavior. For this reason, most 
materials, which are not regular insulators, are metals. Those rare 
materials which have values of Na appreciably greater than 0.3 (n > 2.2) 
owe their existence to a rather delicate balance of forces. 
The third-order energy is 

w 3 = (1| k |1> = -(^)\(E-fi) 2 ) 2 e\{E-nr . (28) 

In most cases a is very nearly isotropic and we have 

** = -«* - (H)« £ - w (2!,) 



so that 



h = -—{{E-ftf). (30) 

na„e 
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With N molecules in unit volume this gives a nonlinear coefficient 

a u c n a e n 

where 

T iik = <£,£,&> = (RtRfa) - (R<)(RM - (RMRJii) 

- QttfPA) + 2(R i )(R i )(R k ). (32) 

Equation (31) is the central result of this section. It expresses d iik in 
terms of the linear (corrected for local fields) susceptibility x and a 
cubic moment (third-order semi-invariant) of the electronic distribution 
in the ground state. 

If we neglect overlap and, for simplicity, also assume that the electrons 
are uncorrelated so that T iik = nt iik where t iik refers to a single electron 
we have 

*** = ^ Un (33) 

and Tat is now apart from numerical factors the octupole moment of 
the charge distribution. If the electron density in the molecule is p(r) 

T iik = JJJ WMz) d 3 r. (34) 

If we account for local fields through a Lorentz correction the correct 
value of x to insert in (33) is obtained from 

n 2 - 1 4t 

^+2 = ¥ X (35) 

and the observed value of d ijk (see Section V) is 

d iik = y — g — ) d iik . (30) 

At first sight (33) seems to imply that d is proportional to x in conflict 
with Miller's rule. However, x depends on N/n(ft 2 ) 2 tt nN{f 2 ) 2 and N is 
inversely proportional to (r 2 )* so that x ^ nr, while d ^ r*. Thus, d is 
in fact more nearly proportional to x 4 than x- 

We now consider as an example, the tetrahedral molecule methane 
CH4 , which crystallizes in the tetrahedral space group F43wi with a 
lattice parameter 5t6 A and a molar volume of 32 ccs. If we take Car- 
tesian axes along the sides of the cubic cell, the bonds point in the 111, 
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and tetrahedrally related, ill, 111, 111, directions. From symmetry, 
there is only one independent component of d ijk , in which all the sub- 
scripts are unequal. 

The shortest c-c distance is 4.2 A and from the size of the free molecule 
we conclude that overlap is unimportant. 

Turner, Saturno, Hauk and Parr 10 have used one center wave functions 
to calculate the electronic density p in the molecule. From this we can 
obtain t 123 using (34). 

The result is 

/,,., - 0.5 X 10" 24 cm 3 (37) 

and this is not very sensitive to the limits of integration. The experi- 
mental molar susceptibility of CH 4 is 1.6 ccs and so x = 0-05. Thus, if 
we neglect correlations between the eight valence electrons we have 
from (33) 

d 123 = 3 X 10" 9 esu. (38) 

In a similar way, neglecting correlations, we can calculate the molar 
susceptibility from (27). Turner et al's charge density leads to 

(x 2 ) = 3.3 X 10" 17 cm 2 

and so, with eight valence electrons, we obtain a = 6.5 X 10" and 
X = a molar susceptibility of 3.9 ccs, rather over twice the experimental 

value. 

This is a clear indication that the electrons are correlated. However, 
the correlation enters twice in x but only once in d (since we have ex- 
pressed d in terms of the experimental x)- Thus, d problaby lies between 
2 X 10" 9 and 3 X 10" 9 esu. 

To see whether 3 X 10"° esu is a reasonable value for d 123 we compute 
the Miller reduced tensor d/x = A I23 = 24 X 10" 6 esu. 

This is quite exceptionally high. Most materials have allowed com- 
ponents of A near 3 X 10"° esu and only one coefficient in LiNb0 3 
(9 X 10" 6 ) and A 123 in hexamethylene tetramine (15 X 10" 6 ) approach 
this value. 

However, we believe it is in fact not far wrong. In most materials 
geometric factors conspire to reduce d by various factors of cos and 
the atomic groups are in the first instance less aspherical than CH 4 . 
In CH 4 the effects of every electron are directly additive. 

Hexamethylene tetramine (HMT), the other exception to Miller's 
rule, is, like methane, a tetrahedral molecule N 4 (CH_ 2 ) 6 _in a_tetrahedral 
/43m crystal. The 4 nitrogen atoms form the 111 111, 111, 111, corners 



NONLINEAR OPTICAL COEFFICIENTS 921 

of a regular tetrahedron and the CH 2 groups occupy the edges but the 
N-C-N bonds are bent outwards in such a way that all the angles are 
very nearly tetrahedral. The carbon atoms occupy the six sites 2, 0, 
etc. (see Kitaigorodskii 11 ). 

The refractivity of the molecule as a whole can be very satisfactorily 
accounted for by a system of additive bond refractions. (See LeFevre 12 
for a review of bond refractions.) The three basic units are 12 C-H 
bonds pointing in the T, T, T, and related directions, 4 nonbonding orbitals 
on the nitrogen atoms pointing in the 111 and related directions, and 
12 N-C bonds in the 111 directions. 

Since the refractivities are additive, these units appear to act inde- 
pendently in determining the molar refractivity. Le Fevre (loc. cit.) gives 
values of R = (Air/3)La (where L is Avagadro's number) of 2.8 ccs for 
each unbonded nitrogen pair, 1.65 ccs for each C-H bond and 0.62 ccs 
for each N-C bond. Thus, the N-C bonds make a rather small con- 
tribution to x, and probably even less to d since they have an approxi- 
mate inversion centre at the centre of the bond (C and N are similar 
atoms as compared with C and H). We therefore neglect them. 

The 12 CH bonds in the ITT direction are roughly equivalent to 3 
methane molecules in the molecular volume 105 ccs, and further the 
electrons will be less correlated than in methane. Thus, their contribution 
to d 123 is 

rffll = - t!t-¥-3 X 10 -9 = -2.7 X 10~ 9 esu. (39) 

To calculate the effect of the nonbonding nitrogen electrons we assume 
that they occupy SP 3 hybrid orbitals directed along 111 etc. with Slater 
radial wave functions Ar exp ( — 2.5r/2o„). It is then straightforward 
to show that for one electron 

t*» = -0.055 X 10~ 24 cm 3 . (40) 

The contribution of the 4 nitrogen atoms to x is 

N 3 4 X 2.8 nn9 „ 
X = 47^05- = °-° 255 

and so 

<& = -1.7 X HT'esu. (41) 

Thus, the total value of d 123 is —4.5 X 10" 9 esu. This could be slightly 
increased by the effects of atomic overlap, and possibly by contributions 
from the N-C bonds. It could be either increased or decreased by elec- 
tron correlations on individual CH 2 groups. The experimental values 
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for the electro-optic effect Heilmeyer" and second harmonic generation, 
Heilmeyer, Ockman, Braunstein, and Kramer, 14 when corrected for local 
field, effects using a Lorentz factor, both give d = ±8.2 X 10~ 9 esu. 
Thus, our calculation is within a factor 2 of the observed value. 

This method is therefore capable, in simple cases, of predicting the 
magnitude of d rather successfully. Moreover, the experimental value 
of d for HMT suggests that we were correct in assuming that CH 4 
will have an anomalously large reduced tensor A 123 • 

The fact that the division of a complex molecule such as HMT into 
simple components leads to a reasonable value for d leads us to hope 
that a similar procedure will be possible in other cases. It might then 
be possible to assign empirical values of d to basic components such as 
the C-H bond or the N: nonbonding pair, and to combine these ad- 
ditively (with a proper attention to geometry) to predict the values of 
d for even more complex molecules. This would not be surprising since 
a similar procedure (see LeFevre loc. cit.) works very satisfactorily for 
the linear susceptibilities. 

It is then obvious that large nonlinear effects will only result, if the 
molecule contains polarizable groups disposed in an arrangement which 
results in a ground state of, far from inversion, symmetry. The large 
value of A in HMT results from the fortunate coincidence that the 
most polarizable components are themselves strongly asymmetric and 
so oriented that their effects are additive. The much smaller values of 
A commonly observed can then be explained as partly due to no group 
in the crystal being quite so asymmetric as N: or CH 2 in HMT and 
partly due to unfavorable geometric relations between the groups. For 
example, if our approach is correct we should expect the analogous 
compound adamantane (CH) 4 (CH 2 ) 6 in which the nitrogens are re- 
placed by CH groups with the CH bond along 111, etc. to have a d 123 
appropriate to 2 ( = 3 — 1) CH 4 molecules in 105 ccs, i.e., d i23 ~ 2 X 10~ 9 
esu or about half the value for HMT. 

Exceptionally small values of A will occur in materials where most 
of the molecule possesses local inversion symmetry, so that only a 
fraction of the molecule contributes to d, while the whole molecule 
contributes to x- We shall consider an example of this in a later section. 

Overlap between adjacent molecules is neccessarily bound to lend 
further uncertainty to the calculation in materials with a pronounced 
band structure, but it seems possible that rough approximations should 
be obtainable from, for example, the relation between bandgap and 
the corresponding separation in the isolated atoms. In fact, since what 
we actually require is T iik /n, which, if the electrons are uncorrelated, 
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is simply 

/Cell 



/Cell 



we may expect that this factor will, to some extent, be self-cancelling. 
Finally, we may remark that very much better approximations to 
d iik can obviously be made if we know the ground state wave function 
explicitly and also use more sophisticated trial wave functions in the 
variational calculation. It is at first sight surprising that a knowledge 
of the ground state wave function alone is sufficient to determine x and 
d which, in the more usual treatments involve the properties of excited 
states. However, we should remember that a knowledge of the exact 
ground state wave function is, except in pathological circumstances, 
sufficient to determine the unperturbed Hamiltonian; thus, the whole 
spectrum of states. 

III. THE CLASSICAL ANHARMONIC OSCILLATOR MODEL 

Although the considerations of the preceding section are sufficient 
to determine the magnitude of d at low frequencies, they offer little 
guide to the variation of d with frequency and, if recast in terms of time 
dependent perturbation theory they tend to lose their attractive sim- 
plicity. In the next section we shall show that a more familiar form of 
time dependent theory leads to results which can be represented in 
terms of a classical anharmonic oscillator model. Here, we discuss the 
properties of the model itself. 

We assume that unit volume of the material contains N optical elec- 
trons which move in a potential 

V = imtfa* + Vntf&fZu , (43) 

where a sum over repeated subscripts is implied. The potential V iik 
obviously satisfies V iik = V iki , etc. 

In a field E^e 0t the equation of motion is 



m m 

and the linear response obtained by neglecting V iik is 



(44) 



-E% i0 ' 
tf ' - S^- (45) 
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There will be a similar response to a field E\e iyt and, if we introduce 
these responses back into the nonlinear term in (44) we obtain a re- 
sponse at the sum frequency a = |9 + 7 given by 

x <« = _3£<^! 1_ 1^ 1^ (E »E1 + ElBfa'" . (46) 
m fi, — a u,- — p il k — 7 

The resulting polarization is Nex\ 2) and so the nonlinear coefficient is 
, 7 ZV iik Ne* 1 L__ I . (47) 

ttk = — 5 ~2 2 ,-.2 a 2 /-k2 2 V**> 



<ft 



tf - a 2 0? - ($ 2 02 - 7 s 



Thus, the symmetry of d iih mimics that of V iih if we neglect the res- 
onance denominators. 

The linear susceptibility obtained from (45) is the familiar expression 

x ...^ * (48) 

m Qi — u) 

and so if we express d°f k as 

ffff = Xli&xl^m (49) 

the reduced Miller tensor is 

a = JL&tt (50) 

A *'* iVV ' 

which is frequency independent and has the same symmetry as V ijk . 
If we assume that V iik is electrostatic in origin its order of magnitude 
will be e 2 /d 4 where d is an atomic spacing and we shall also have Nd 3 « 1. 
Thus, 

I A,,, I « 8 f (5D 

With d equal to 2 A this is 2.5 X 10" 6 esu, about the mean value of A 
for most materials. In a later section, we shall give another estimate of A. 

The potential V iik XiXjX k distorts the shape of the ground state of 
the harmonic oscillator and as a result the system acquires a cubic 
moment t iik [defined in (32)] which we now calculate. 

Let I 0) represent the unperturbed ground state wave function in 
the absence of the anharmonic term and | p) be an excited state, then 
the perturbed wave function is 

l> = |o>-£'^pl P >. (52) 
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The expectation values of even operators such as (x*), (z.z,) are un- 
altered by V, while the expectation value of an odd operator such as 
x, or XiXjX k is given by 

<w.) = -2 E' <°l '**. WpI V |0> . 

It will suffice if we calculate t iik with i 9^ j ^ k. Since (x<rc,-) = if 
i 7* j we only require (x { X;X k ) and contributions to this come only from 
the 6 = 3! terms in V with i ^ j *£ k. The only state which contributes 
to the sum is | p) = | 1, 1, 1) with an energy fc(G, + fl 2 + 0,). The matrix 
element is 

<o| x,«, |in> = (A) ! (_I_y 

and so 

*» = <***> = -12(55] Ao,o,a s (o 1 + o I + n 3 )' 

It is straightforward to show that a similar result 

/ = -2 (h-Y Yjh nn 

" k 2 W KQiQiQk(Qi + 0* + «*) ^ 

holds for all the components of t ilk . 

If we substitute this relation in (47) and take the limit as a(3y — » 
we obtain 

diik = 2 ^e tiik • (55) 

This is twice the value obtained in (33) because there we treated t iik 
as a fixed property of the ground state which was then perturbed by E; 
whereas here we have considered an even ground state perturbed by E 
and a fixed potential. 

Thus, if the real system has a cubic moment t ijk in the ground state, 
the equivalent anharmonic oscillator model requires an anharmonic 
potential 



V' iik = -^S t iik (56) 



and this will result in a cubic moment t' iik = \t iik in the oscillator ground 
state. 

In the real crystal t ijk may be an accessible quantity. It obviously is 
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in molecular crystals of strongly covalent compounds such as CH 4 . 
But, in ionic crystals it may be more sensible to consider the ions as 
spheres perturbed by a crystal potential V' iik . In a later section we 
shall see that there is a simple relation between the model potential 

and V e m . 

The classical anharmonic oscillator model has previously been used 
by Bloembergen, 15 Garrett and Robinson 16 and Kurtz 17 to give a qualita- 
tive account of nonlinear phenomena. The latter authors also discuss 
in some detail its relation to Miller's rule. 

Obviously, the model is the nonlinear analogue of the classical har- 
monic oscillator model used with such success for the last 60 years in 
the discussion of linear behavior such as dispersion, and, just as the 
harmonic model is directly related to the results of a quantum mechan- 
ical treatment, we may expect the anharmonic oscillator to have a 
similar basis. In the next section we explore this relation. 

IV. TIME DEPENDENT QUANTAL TREATMENT 

A number of authors Bloembergen, 15 Armstrong, Bloembergen, 
Ducuing and Pershan, 18 Butcher and McLean, 19 Kelley, 20 Cheng and 
Miller 21 and Ward 22 have given rigorous quantal treatments of optical 
nonlinearities in solids. We select an expression due to Armstrong, 
et al (loc. cit.) which expresses the nonlinear coefficients in terms of the 
energies hw v of excited states and the matrix elements (0 | z, | p), 
(V I %> I Q.)} etc - of tne dipole operator between states. The ground state 
is <0 |. 

This expression is valid, either for an assembly of N isolated atoms 
in unit volume or, in the dipole approximation, for a real solid where 
the wave functions overlap. In the latter case, the solid must be divided 
into cells of the periodic lattice, and N is then the density of cells, while 
the matrix elements are to be evaluated only over the interior of a cell. 
The periodicity of the lattice ensures that contributions from parts of 
the wave function outside a cell cancel in the crystal as a whole. 

To avoid a plethora of subscripts we let each of x, y, and z serve to 
represent any one of the components and we can then write the expres- 
sion for d as 

iaBy Ne 3 -^ x^ / ggga + "T 

d - - T ? ? !"■**• (4 -T)« - J) 

co p co q + /3a . (Jpttf, — 7 I , rr , 

+ iwfcA. (w ; _ ,*■)«- a 8 ) + ZMJvo <& - y)(4 - f)) { V 
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This expression vanishes if the states | 0), | p), etc. have a definite parity, 
its value therefore depends on the existence of matrix elements whose 
presence is contingent on the absence of inversion symmetry. For this 
reason, it is almost impossible to make an informed guess about its 
magnitude or behavior. 
An analogous expression for the linear susceptibility is 

_ 2A r e 2 y, w p x np y pn 
h „ u p — io 

and in both expressions an operator x is to be understood as the total 
operator for the contents of a cell, i.e., the sum of the individual electron 
operators. Of course, we can neglect the core (nonvalence) electrons on 
the grounds that they are too tightly bound to contribute to the optical 
properties. 

A familiar approximation to x is obtained if we note that in (58) the 
variation of the summand with | p) is almost exclusively due to the 
matrix elements. These not only obey selection rules, but also decrease 
rapidly in magnitude as the state | p) increases in energy, and therefore 
overlaps the ground state less and less. For example, in the H atom with 
a IS ground state the matrix element x op vanishes unless p is one of 
the states 2P, 3P, etc. Moreover, as we go from the 2P state to the SP 
state x op x po decreases by over a hundredfold. At the same time, a> p 
changes by less than 30 percent. Thus, except near a resonance, we can 
treat o>„ as a constant ft, somewhere near the first allowed transition and 
write (3.2) as 

X " = h(tf - J) ^ XovVv ° (59) 

where the primed sum excludes p = 0. Now 

£' XopVpo = 2 Xo P y PO - x 00 y 00 = (xy) 00 - x 00 y o „ 

- <(* - V»(y - (y))) (60) 

where a ( ) denotes a ground state expectation value. Thus, 

2A t p 2 Q 
& = , |(tf _ nsj <(* - <*»(?/ - (2/»>. (01) 

We shall not pursue the further manipulations of (61) using the sum 
rule which lead back to (27) but we remark that in many cases a form 
such as (61) for x, involving a single Sellmcier or classical oscillator 
term, gives an excellent account of optical dispersion, and that when 
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applied to the hydrogen atom, with hQ set equal to |e 2 /«» the IS — 2P 
energy, it leads to a value of x at low frequencies 

1 6 3 

X = irflo 

which exceeds the correct value 4.5 a 3 by 32/27 or 18 percent. 

Before we can adopt a similar procedure with the nonlinear coefficient 
we must first satisfy ourselves that there is no essential difference be- 
tween a sum with three matrix elements and one with two. In x all 
matrix elements terminate on | 0) but in d it is quite possible in a term 
such as x ep y vll z Q , , with p & q corresponding to highly excited states, of 
great spatial extent, that the term y PQ may be large enough to compensate 
for the smallness of x op z QO . If this were the case it would be possible for 
the exact value of the sum to depend critically on cancellations between 
large terms involving highly excited states. The removal of the frequen- 
cies w p , etc. as a single average would then have disastrous results on 
the sum. 

We will advance three arguments why this is unlikely. Consider first 
an even higher-order calculation, that of the fourth-order Stark shift 
of the ground state of atomic hydrogen due to a field F. In atomic 
units this is given by an exact calculation (Dalgarno 23 ) as 

w w = -m^F* « -5GF\ (62) 

We can also express W (4) (Dalgarno, loc. cit.) as 

W (4) _ V"V V"»/ y»' X BV X PQ X QT X rn xp, X op X v „ y>, X np X p „ _ 

w ~ r r r cwr + r «. ^ <£ 

Our procedure treats o) p u q and w T as a single constant ft and leads to 
W w = -| 4 5 {((x - <x» 4 > - 2({x - (x)) 7 ) 2 ). (63) 

For the H atom (x) = 0, <x 2 > - 1 a.u. and (x*) = 9/2 a.u. so that, if 
we set ft = 3/8 a.u., the IS - 2P energy difference 

pp-<« = __iffopi ~ _48f 4 . (64) 

This is close to the correct result (36), but despite the fact that we have 
taken the lowest possible value of ft it is too small. This is a clear indica- 
tion that some cancellation of higher terms, which we have aggravated 
by our cavalier treatment of co p , etc., is occurring. This is not surprising 
for, if in the triple sum we consider the lowest possible sequence of levels 
IS 2P 2S 2P IS for which o>„ = a> a = co r = ft the product of the matrix 
elements is 5 a.u. while for the sequence IS SP SS 8P IS the product 
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is 10 a.u. made up of a contribution of 0.0033 from the two IS 8P 
elements and 3000 from the 8S 8P elements. 

However, this is not quite so serious as it appears, for in a real solid 
no matrix element can exceed the linear dimensions of a cell say 5 a.u. 
so that the product in the low transition would remain at 5 a.u. while 
the product for the upper transition would be reduced to 0.08. 

Our final argument is empirical. If cancellations between large terms 
are critically important, the relevant feature of our procedure is the 
change in the ratio a> p /a> g it causes for highly excited neighboring states. 
In hydrogen the ratio of the IS — SP energy to the IS — IP energy 
is 1.005 and we replace this by unity. In a time dependent theory res- 
onance denominators appear, and if the sum is really so critically 
balanced, we expect the observed quantity, in this case the hyper- 
polarizability, to vary rapidly with frequency when w 2 za 0.005 fl 2 , 
i.e., at a frequency 10 times lower than the first absorption edge. In 
nonlinear optics, no such variation is observed until one of the frequen- 
cies approaches much more closely (about 70 percent) to the absorption 
edge (Chang, Duelling, and Bloembergen 24 ). 

Taken together these arguments give us reasonable grounds for 
hoping that the sums will not bite us if we remove a> p , etc. from under 
the summation sign. 

In the sum in (57) there is no restriction on p or q, in particular terms 
with either p = Oor q = occur. These will lead to trouble if we attempt 
to approximate the sums as they stand. We therefore first segregate 
all such terms. Let { } denote the entire summand in (57), then 

p 4 p « a r w r — 7 7 r w r — a 

_ a y Z er X„ fiy, y or z ro 

■ £ ^ X or y ro . 7 VV Z -r X ro /r -\ 

+ Zoo - z, -2 ^ + y no o 2^ -? — -i- (6o) 

7 r w r — p p r OJ r — 7 

Three single sums remain unprimed, but because a = + 7 the terms 
with r = cancel and so we may regard all the sums as primed. 

We can now remove w p , a>„ and w r as a single average £2, and this leads 
to an expression containing terms such as 

V <i pa r 

- Z O0 2 XorVro + 2x„„1/ uo Z„ n . 
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Each of the sums on the right is now a ground state expectation value. 
When all the terms are collected together we obtain 

d -n _ Ml «W + fr - «') Hxuz) _ (x)(yz) _ mzx) 



- (z)(xy) + 2(.r)(//)( 2 )| (66) 



which we can also write as 



<%!: - ^D(Q,a,fi f y)T w (67) 

in terms of the, by now, familiar cubic moment T IV2 . This expression 
bears an obvious resemblance to (61) for x- 

Our expression (66) or (67) would be very nearly exact if all the optical 
levels had very nearly the same energy. It would then correspond to 
the fictitious two level system (see Refs. 15, 16, 18) often used to oblit- 
erate some of the intractable features of (57). Unlike this model, how- 
ever, our expression retains the geometry of the system implicit in the 
selection and sum rules. 

Equation (66) is possibly valid up to a frequency where one of a, /3, 
or 7 approaches the first allowed transition frequency. At somewhat 
lower frequencies, it is legitimate to drop the term 0y — a in the numer- 
ator. This then allows us to make a further generalization at no increase 
in complexity. 

By removing o>„ and u) Q from (57) as a single average we have tacitly 
neglected the possibility that the system might be birefringent. We 
can remedy this by noting that in (57) each frequency u p or «, 
is uniquely associated with a matrix element such as x op or z Q0 which 
terminates on the ground state | 0) and therefore also appears in x- 
Thus, we can consistently introduce three averages Q, x , S2„ , and ft, as- 
sociated with correspondingly polarized transitions. If we follow this 
process through all its tedious ramifications, we find that, except in 
the term /3? — a which we are omitting, it leads to the surprisingly 
simple result that D(Sl,a,f3,y) is replaced by 

nm a \ MMQ, + Q, + gj , m 

m «. 0. y) = (fl j _ ;-) ( q{ _ rtoaZ - 7) m 



Thus, 



lol3y _ Ne A UMUn, + 0, + fl,) r m) 

d «" - tf (tf - a 2 )(nl - ?){tf. - r) J ' y - ' K ' 
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If we compare this with the result for the classical anharmonic oscillator 
obtained by combining (54) with (47) we see that they are identical 
except for a factor 2 which once again arises because in one case we 
assumed that T xyi was a fixed parameter while in the other it was V xyt . 

We now see that the classical model is equivalent to the quantal 
treatment, except near a resonance, in the following sense. 

If we construct the model, by choosing Q x , Q y , and Q, to give the 
correct linear properties then we must choose the anharmonic term in 
the potential to produce a cubic moment in the ground state of the model 
equal to £ the corresponding moment in the real system. The dynamical 
properties of the two systems are then equivalent and the model can 
be used to treat more complicated systems where the quantal treatment 
is too difficult. 

We now consider the relation of V xul to the actual potential respon- 
sible for the existence of T xvt . Obviously, the relation is obtained by 
requiring that both potentials yield the same cubic moment, one in the 
model, the other in the real system. In this case, there will be no factor 
of 2. 

For simplicity, we consider only a system which is isotropic before 
the application of the anharmonic potential. Further, we restrict our- 
selves to atoms in which there is only one valence electron. Our results 
will, however, be directly applicable to atoms with more electrons if 
we can neglect electron correlations. 

We already have an expression for the oscillator (54) which we can 
write as 

r...--*t%£. (70) 

where a>„ , the classical frequency, also corresponds to the first allowed 
transition, and 

" - fey < 7 » 

is a measure of the extent of the system in one dimension. The direct 
proportionality between T ijk and the corresponding component of 
V iik occurs because the oscillator Schrodinger equation separates in 
Cartesian coordinates. In general, as we shall show, it will only hold if 
V = V iik XiXjX k , the crystal potential, satisfies Laplace's equation. 
We will consider a more general potential of the form 

V = Z V:^P m t (e, v ), (72) 
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where P" is an associated Legendre polynomial normalized to unity. 
This potential satisfies Laplace's equation only if n = I. 

If the unperturbed ground state wave function is \p the first-order 
correction ^ due to V satisfies 

(//„ --£?„)*, + (V - A',)*,, - 0. 

Since T is an odd moment we need only consider odd terms in V (in 
fact only / = 1 and / = 3) and for these E x vanishes since f has definite 
parity. 
We let 

fa = ty (73) 

and then 

(H Q - E )f+ Q = - Vf 

but, since 






this leads to 



V7 + 2V/-V log* = y~ Y. (74) 

Now ^o is a function of r alone and so we can write 

/ = Z V; ianl (r)P1(d,ip), (75) 

where a n ,(r), which does not depend on m, satisfies 

1 * , da _ 1(1 + 1) , „ ?a fl leg +„ = 2iv 

? dr F dr ' r 2 " + " dr dr ' tr ' 

The perturbed ground state is therefore, 

* = {1 + E V:,Qn l (r)F m i(0,v))MT) (77) 

ii I m 

and in this new ground state we can easily evaluate expectation values 
such as 

<r'JV> = £ FSArt , (78) 



where 



ft*- f° r'a nX (r)ti(ry dr. (79) 

Jo 
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In evaluating T iik we shall need (a:,), (x*), and (XiX t x k ). The even 
moments are unchanged by V and we obtain the odd moments by ex- 
panding x { and XiXjX k in terms of Legendre polynomials and powers 
of r. 

We omit most of the gruesome details of the calculation, and further 
restrict V to contain only terms of the type 

v = v iik x iXi x k + x^ = £ f 3 vp : 7 + vryp? + m**? . (so) 

The term in F 3 ™ which does not satisfy Laplace's equation is necessary 
to obtain the most general form of the cubic part of the potential 
VijkXiXjXi . This contains 10 independent parameters while P% has only 
7. The missing 3 are supplied by P'". 
If this term is absent, we have 

V 2 F = (81) 

and then 

S t = V iU + Via + V ikk = 0. (82) 

With all the terms present we obtain 

<:<-,) = fanSi + tfinXi (83) 

<a3) = JtfmVm + 3{t&/3 83 . - rh&^Si + ifi ni X t (84) 

(*,*•) = Aft M 7„, + |M, 31 - Th0 3aa )S t + M,,!,. (8f>) 

<.r,.r,,r t .) - Aj8 w V, /4 , (86) 

where it is to be understood that i '^ j ^ k. 
If we let 

t = (•<•;), (87) 

then since <.r,.c,) = 0, i i* j, and (x t )(x t )(x k ) is third order in V, we 
obtain 

7' fll = Afta.F,„ + 3 {Aft,,, - rririS 3: , a - §7/3,3.15* 

+ 3{A/3 3I1 - fr/3 m }X, (88) 

7\„ = t^j8 333 F„, 4- |A/3 331 - r^i8 3 33 - hAu}ft 

+ {Aftn ~ MiiiW (89) 
7',,, = Aft„7i, 4 . (90) 
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For a harmonic oscillator 



a - — q- — 
/ * 333 ~ d ° tua a ' 



^ = ~ S5 t' ^ = - U t 



, - _a_ ,- a 



7 = a 



A,, = -3 



hu„ 



(91) 



and it is easy to check that the coefficients of 5, and X { vanish, so that 
we recover (70). 

If V satisfies Laplace's equation S { = and in the absence of an 
internal field X t every component is given by 

T iik = AA»F W . (92) 

Thus, in this case T ijk and the reduced Miller tensor A,,* have the same 
symmetry as V iik . Therefore, since <$»,• = Owe have 

A,.-,- + A,-,-,- + A,„ = 0. (93) 

If i is an axis of 3-fold or higher symmetry, A,,, = A fU and so, for ex- 
ample, 

A 33 3 = -2A 311 . (94) 

This relation is rather well obeyed by the coefficients for the 6-mm 
crystals listed in the Table I. Signs are available only for the electro-optic 



Table I 



Material Wavelength m Ajja X 10 6 esu Jim X 10' Mil 



Ratio 



Linear Electro-optic 



ZnO 

ZnS 
(MS 



optical 
optical 
optical 



1.6 
O.fl 
1.2 



-0.8 
- . 45 
-0.55 



Second Harmonic 



-2.1 
-2.0 

-2.2 



ZnO 


LOG 


3. 3 


1.1 


±3.0 


ZnS 


10.6 


4.9 


2.45 


±2.0 


CdS 


1 . 06 


3.2 


1.6 


±2.0 


CdS 


10.6 


5.4 


3.3 


±1.6 


CdSe 


10.6 


4.8 


2.4 


±2.0 
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coefficients and so the s.h.g. results represent moduli only. References 
to the experimental data are given in conjunction with later tables. 
In the case of the electro-optic data, the experimental figure is for A 113 
and we have assumed that Kleinman's rule (Kleinman") holds and that, 
this is equal to A 31 , . Except for s.h.g. in ZnO the ratio is —2 within 
the experimental error. 

If, on the other hand, the sole perturbation in V is the field .Y, , we 
have 

Z\„ - ZT W = §!&,, - 57/3,n) (95) 

and the expected ratio is +3. In crystals where both terms occur in 
V with arbitrary strength, any value of the ratio is possible. This is 
observed in the ferro-electric crystals BaTi0 3 ratio +£ and LiNbO, 
where it is +1.7 for the electro-optic effect and ±11 for s.h.g. It is 
perhaps somewhat surprising that the ratio is so exactly —2 in the 
6-mm crystals since this is a polar point group and an internal field 
X a is not forbidden by symmetry. 

If V does not satisfy Laplace's equation, (it need only satisfy Poisson's 
equation) there is no direct relation between the components of T iik 
and those of V iik even in the absence of a field X, , although, since 
xyz is a spherical harmonic, we still have 

7'm = WjSssa ^123 • (9f>) 

However, since the coefficients of 8 t vanish for the harmonic oscillator 
we may expect them to be small in other cases. We gain some support 
for this view by considering the hydrogen atom for which 



105V a„ a 35385 al 1485 at 

Pa 1 1 — 



, (97) 



8 / foo, ' ™' 256 hw t ' "'" 04 /to, 

_1305^ _ _315^ 81 a! 

P,M "" 128 fia>, ' Pi * ,B " l(i k, ' 0, "~~32^ 1 

y = o'i 

where as usual a = ti'/mc' and Jta, = 3/8(e7a<.) i« the first allowed 
transition (IS — 2P) energy. 

The coefficient of V iik in each term of T iik is then - 315/1 Q,{al/tiu} x ) 
while the coefficient of £, in T iU is a factor 23/200 smaller. In T ifi it 
is 23/()00 smaller. Thus, in the absence of X, the non-Laplacean terms 
in 1' cause no more than an 11 percent departure from the relation 

r "' -_ i5"c: F "*- (9 « 
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Since we expect the dominant terms in V to satisfy Laplace's equation 
it appears that T iik , A,-,* and the model potential V' ijk will be very 
nearly proportional to the corresponding terms in V. 

The potential V required in the model is related to the crystal 
potential by 

For a hydrogen atom this gives V' iik * 5V iik thus, insofar as real atoms 
behave like hydrogen atoms, a model with the same spatial extent 
a « a and the same first allowed transition o> « wi will require a 
potential roughly five times as strong as the actual potential. This 
reflects the obvious fact that a harmonic oscillator is a stiffer system 
with more sharply localized (\f/ « e~ r ') wave functions than an atom 
(* « e' r ). 

We have now shown that, with an appropriate choice of parameters 
a classical anharmonic oscillator model is a very good approximation 
to the intrinsic electronic nonlinearities of real systems. 

In the next section, we use the model to consider the effect of lattice 
polarizability which we have so far neglected. 

V. LOCAL FIELDS AND LATTICE POLARIZATION 

We have already remarked in the introduction that the seat of the 
nonlinearities resides in the electronic motion. It is, however, consider- 
ably modified by local field corrections and in the case of optical rectifica- 
tion and the linear electro-optic effect by lattice polarization. 

Miller's rule states that dgf is proportional to the product of the 
observed linear susceptibilities X," , etc. at the appropriate frequencies. 
If one of these is dc we are to take the actual dc susceptibility and not 
the extrapolated long wavelength limit of the optical susceptibility. 

At first sight, it seems plausible that this is simply the effect of in- 
ternal fields, which cause the local field experienced by an atom to be 
greater than the applied field. We now examine this hypothesis and 
show that it is inadequate. 

Microscopic calculations yield the polarization of single atoms due to 
local fields. In the linear case, if we have N atoms per unit volume of 
polarizability a 

P = Ned*], (100) 

and the local field is related to the applied field E by 

Et= E + TP. (101) 
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In some cases the Lorentz value of r = -iw/3 is applicable and we then 
obtain the well-known relation between the refractive index n, or the 
dielectric constant e = n 2 and a. 

n 2 - 1 4ir ., It , . 

where V is the molar volume and R is the molar refractivity. 
In general, 

P - Na(E + TP) = t *° Na B (103) 

and the observed susceptibility is 

while 

E L = (1 + V X )E = YZ^gi' (105) 

In nonlinear optics the two driving fields E* and E\ are obviously 
modified according to (105) but, as Armstrong, Bloembergen, Ducuing 
and Pershan 18 have shown, there is a further factor in P. This arises 
because the nonlinear polarization 

V, = d?A 7 (2$,„..,(^),„..i , (10(5) 

produced directly on the atoms, further polarizes the surrounding 
medium. 



We have 



so that 



P; = p? -f- TNaPt , (107) 



P? = r=TN~a = (I + TX " )P * • (l08) 

Thus, if d'ft is the (calculated) intrinsic coefficient, the observed 
coefficient is 

Dfff = (i + r x ?.)(i + r x ?,)(i + Txi^n??. (109) 

Therefore, even if d does not vary with x, D will do so. This is, however, 
not enough to explain the observed variation of D with x- For example, 
in semiconductors it is very likely that V is small if not zero and yet 
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the measured values of D appear to obey Miller's rule and be propor- 
tional to x 3 - Thus, the intrinsic coefficient d itself must have a similar 
dependence on X, 
If we write 

in terms of the measured susceptibilities (i.e., ri — 1), which is the con- 
tent of Miller's rule, and then use (104) to express D in terms of the 
atomic polarizabilities we obtain 

DJff = (i + r x .")(i + r x f,)(l + TxtotftidALii* ("« 

so that from (109) 

dgf _ Na'uNo^iNalu^u . (112) 

Thus, the reduced tensor is the same whether or not we apply Local 
field corrections as long as we do it consistently. To obtain a more or 
less constant value of A we must have d varying as a. 

Since A for NH 4 H 2 P0 4 derived from the purely optical s.h.g. effect 
agrees with A from the quasi-static electro-optic effect to within 10 per- 
cent, although the values of d differ by a factor of 12 and in BaTi0 3 
the two values of A au are within a factor 2 while the d's differ by 300 
it is clear that lattice polarization has a direct effect in d not described 
by local field terms. 

We repeat that optical nonlinearities have an electronic origin. 
Electrons in atoms do not move in a harmonic potential. Second har- 
monic generation, which can only involve electronic motion, is much 
the same in covalent organic materials, ionic crystals and ferro-electrics. 
Large values of d 1 " are associated exclusively with large refractive in- 
dices. Thus, nonlinearities in the ionic motion play a secondary role 
in nonlinear optics; however important they may be in determining 
the ferro-electric properties. 

We shall attempt to construct a model, just sufficiently general to 
exhibit the gross features of ferro-electric behavior, and show that it 
modifies the nonlinear optical behavior exactly as predicted by Miller's 
rule. The model is not put forward as an explanation of ferro-electricity 
although it has a venerable past in that connection, but as a demonstra- 
tion that a simple system with singular dielectric properties behaves 
in a way consistent with Miller's rule. 

In Fig. 1, we illustrate a moderately realistic one-dimensional model 
in which electrons of mass m are coupled to ions of mass M in a lattice. 
Forces act between like and unlike particles and of these by far the 
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M K MM M K MM M 

Fig. 1 — A realistic nne-ilimcnsioiml model. 

strongest is K mM which is responsible for the electronic optical spectrum. 
The remaining forces determine the lattice spectrum. The important 
nonlinearities are associated with K Mm . The linear behavior of this 
model is formidably complicated and we therefore assume that its 
salient features are already evident in the much simpler model of Fig. 2. 

The electron of mass w, = m is coupled to the ion of mass m 2 = M 
by the force constant k l2 which replaces K mM . It is anharmonic. The 
electron and the ion are also coupled to rigid supports representing the 
rest of the crystal by forces /c, and k 2 . It is as though we had gone di- 
rectly from the Born- Von Kantian theory of specific heats to the Ein- 
stein theory without mentioning Debye. 

Let Xi be the displacement of the electron of charge e, and x 2 that of 
the ion of charge e 2 . We shall assume that the potential energy is 



<P = \hx\ + |fe 2 &2 + fefri - x 2 ) 2 + v 12 (x } - x 2 ) z 



(113) 



so that the anharmonic term is exclusively associated with the "atomic" 
binding of the electron to its parent ion. It will be convenient to define 
v 2 i = — v l2 . The equation of motion in a field E e tpt is then 



0i0t 



nijXi + ki.i'i + k,,(.Vi — x^ + SVi,(Xi — .r,)" = e t E e 



and the linear response neglecting v tj is 

(/,-, - m,/3>, + /,-,.>(<', + e 2 ) 



." ' 



(/,-, + /,-,, - m,fP)(k 2 + k ]2 - m 2 2 ) - kU 
With A T units in unit volume, the polarization is 
P a = Nfrxl" +«»») 



-E*e m . 



(114) 



(115) 
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r//'///////,. 
Fig. 2 — A simplified one-dimensional model. 



and so 



. N e\(k 2 - m,0 2 ) + eljkt - w,/3 2 ) + fc^jgi + c a ) 2 . (1 j fi) 



X (A:, + A- 12 - m l fi 2 )(k 2 + fc ia - m 2 2 ) - fc?„ 

At an optical frequency u well above (kz/mrf the ionic resonance 



. w Uh - 2i (ii7) 

fcl + «12 — Wl,« 

while at flc 

o _ M feggj + feid + ^ i2 ( Ci ~*~ ga) (ii») 

To obtain the sum frequency polarization due to two fields E e ' 
and E y e iyl we substitute the linear responses back into the nonlinear 
term in (114). The result is a nonlinear coefficient 

d*f = -3*W(a)/03)/(7), (119) 

where 

/(a) = ft + t* - m ia 2 )(k 2 + /c 12 - m 2 a>) - kj* ' ^ <""' 
If we express d aflT as x"x X y & we have 

A- -^g(a)g(0)g(y), (121) 

where 
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k a — rriafx — — (A - , — w,a") 

gift) = — . — ^ , etc. (122) 

A L , - m,a 2 + \^j (A-, - m l a z ) + fc I2 (l + fj 

We note first that if, as seems most reasonable, e 2 = — e x then g(a) = 
g(0) = g(y) = 1. In any case at optical frequencies gr(co) « 1 for all 
reasonable values of e 2 /ei and at dc 

lc, - - A:, 

9(0) = - -Ty —, T 2 (123) 

, 2 + (|), + 4 + .) 

which is also near unity if e 2 ^ — e x . Thus, to all intents 

A « -|H (124) 

which is exactly the result obtained by neglecting the ionic motion. 

Thus, A is an intrinsic electronic property and the effect of ionic 
motion is entirely contained in its effect on x- We note, however, that 
in some ferro-electrics, where the departure from inversion symmetry 
is both small and temperature dependent, A will also be temperature 
dependent. 

If /ci , /c 12 , and fc 2 are all positive, the dc susceptibility is greater 
than the low frequency limit of x" but not dramatically so. There is, 
however, no reason why one of these constants should not be negative. 
Negative compliances are familiar in classical mechanics, a well-known 
example is the common automatic door stop which exhibits a positive 
compliance as the door is first opened but a negative compliance when 
the door is almost fully open. The force between atoms as a whole in a 
lattice exhibits a positive compliance but if we separate this force into 
nuclear-nuclear and electron-electron repulsion and nuclear-electron 
attraction, it is quite reasonable to assume that at the equilibrium 
distance the latter component has a negative compliance. 

It is immaterial which term in (113) we take as negative although 
on physical grounds it seems most suitable to take fc, and this is also 
a convenient choice. 

Provided that 

A-,A- L , + fc 12 fc, + k ls k 2 > (125) 
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or 



, aB _ tl fejtil<l (126) 

K2K12 

the system remains in stable equilibrium at x, = x a = 0. 
The natural resonance cox and co 2 of the system satisfy 

mitn&lwl = kik 2 -+- kik 12 + k 2 k yl (127) 

and so as rj — > 1 one of these frequencies — >0. At the same time the dc 
susceptibility (for simplicity we take — e 2 = e Y = e) 

i _ k\ 2 

o Ne 2 k v , + k 2 (ioq\ 

A, 2 i — v 

becomes infinite, while the low-frequency limit of the optical suscepti- 
bility remains finite. 

If ij exceeds unity there is a spontaneous polarization limited only 
by terms such as ux 2 which we have failed to include in <p. 

All this is reminiscent of ferro-electric behavior if 77 is temperature 
dependent and the Curie point corresponds to t\ = 1. 

The inclusion of a term coxt in <p will, in fact, make y temperature 
dependent, for the effect of this term is to replace k 2 by an effective value 
for low-frequency displacements 

k' a w k 2 + 6co^ = A- 2 (l + XT), (129) 

where x\ is the mean square thermal displacement. As a result if r\ n is 
the value at T = we have 



and so 



k 2 ■ -.1 jjjia m m 

_ JVe 2 fcig \fc12 + fc 2 

x := k 12 " x(r- t ) 



(131) 



if we define 2',, as the temperature at which y = 1. This is of course a 
crude approximation to a Curie-Weiss Law. 

By ascribing all the temperature dependence to changes in k 2 , it 
is obvious from (117) that X" is temperature independent. For 17 to 
be equal to unity, it is not necessary for — k x to be of the same magnitude 
as k 12 , all we require [see (126)] is that -fc, be near k 2 . Thus, from (117), 



NONLINEAR OPTICAL COEFFICIENTS 943 

we do not expect any very anomalous values of x" in ferro-electrics, 
except in so far as materials with a high electronic polarizability are 
more likely to be ferro-electric. 

We have now shown that it is possible to incorporate in our model 
features which lead to quite different behavior for the optical and dc 
dielectric constants without either invalidating Miller's rule or even 
changing the value of A which is essentially a purely electronic property. 

We should, therefore, expect the temperature variation of D"fJ to 
correspond to that of x,"x/,Xa* and this is well borne out by the measure- 
ments of Z wicker and Scherrer Jo of the electro-optic coefficients and 
Bass, Franken, and Ward 20 of the optical rectification coefficients in 
the dihydrogen phosphates. Both coefficients are directly proportional 
to the dc dielectric constant which obeys a Curie Weiss Law. 

In KDP there is almost no change of the s.h.g. coefficient (Van de 
Ziel and Bloembergen 27 ) with temperature above or below the Curie 
point, in accord with our expectations, but at the Curie point there is 
a small discontinuous change. In an orthorhombic coordinate system 
rif," and daaj are equal above T„ but below T B , d 2 ," increases and d££ 
decreases while at the same time xn ~~ Xsa decreases and X22 — X33 
increases. With a constant A this is not compatible with Miller's rule. 

However, at the transition there is a change in crystal class in which 
«i increases and a 2 decreases, Jona and Shirane. 28 It is not unreasonable 
to assume that this increases T 3U and decreases T 322 by more than 
enough to compensate for the changes in Xn and X22 • 

vi. miller's rule 

The classical anharmonic oscillator model, which we have shown to 
be a good approximation to the behavior of a real system, leads directly 
to that part of Miller's rule which refers to the geometric properties 
and frequency dependence of the nonlinear coefficients in a single 
material. We have also in (51) advanced a crude argument to show that 
A will not vary much from material to material. 

When we examine the experimental data we shall see that the allowed 
components of A are between 1 X 10" 6 and 6 X 10 -0 esu for most mate- 
rials but that there are a few materials with significantly higher values 
and a number with values as low as 0.1 X 10~° esu. 

In most cases, these exceptional values have a rather simple explana- 
tion and we have therefore to explain a constancy of A to within a 
factor of about 10. 

Neglecting the effects of lattice polarization and local field corrections, 
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which we have shown are irrelevant, the results for the classical an- 
harmonic oscillator model are, from (50) and (56), 

A 3 InjL (132) 

a B ex 

This is also the result from the static perturbation treatment of Section 
II. 

If we use 

x « 4JNT £2- = W G3- , (133) 

a„ a„ 

we arrive at 

_ 243 a, T„, , (134) 

Now the volume occupied by the oscillator is both 1/N and 8 <r 2 )* and 

so 

*m « 10W pyl esu (135) 

where we have inserted numerical values for a and e. This expresses 
A,-,* as the product of a scale factor (r 2 )* and a dimensionless shape 
factor T/r 3 . Whether we assign to each oscillator the volume per valence 
electron, per atom or per group of atoms, (r 2 ) 1 is likely to be between 
0.75 and 3 A; so that A will be sensibly constant near 3 X 10 6 esu, if 
the shape factor is of the order of 0.01 to 0.05. We have from Turner, 
Saturno, Hank and Parr's 10 results for CH 4 a shape factor of 0.05, 
and so this range of shape factors is not unreasonable. It corresponds 
to a linear distortion 0.02 s pa 25 percent. It is also not unreasonable 
that the distortion should be of this general order, wherever it is al- 
lowed by symmetry. We may speculate that much smaller values of T/r* 
would imply very weak interatomic forces and much larger values would 
lead to a structure unstable relative to a more symmetric arrangement. 

Thus, qualitatively, the relative constancy of A reflects relatively 
constant shape factors, although we can hardly claim that this is more 
than a sophisticated form of dimensional analysis. It does, however, 
suggest that A is determined primarily by the geometric properties of 
the molecular and crystal structure. 

Large values of A will occur only when the molecules themselves 
depart considerably from inversion symmetry and are arranged in the 
crystal in such a way that the effects of individual parts of the molecule 
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are additive. Small values of A will occur when sections of the molecule 
have local near inversion symmetry or when their disposition in the 
crystal favors the cancellation of effects from different atomic groupings. 
However the molecules are arranged in the lattice, A will be small if 
the molecules themselves have near inversion symmetry, or consist 
of uncoupled parts with the same property. 

In Tables II, III, and IV, we present 50 values of A derived from 



Table II— Second Harmonic Coefficients 



Units of d 10 _;i esu Units of A 10~ 6 esu 



Material 


Clan 


hx 


di-.s 


| 

A 


Ref. 


HMT = N., (CII 2 )a 


43m 


1.06 


30 


17 










1 


ZnS 


43m 


1.00 


153 


3.5 










■■', 


ZnS 


43m 


10.6 


146 


4.5 










2 


ZnSe 


43m 


1.06 


200 


2.5 










:: 


ZnSe 


43m 


10.6 


370 


6.6 










2 


ZnTe 


43m 


1.06 


660 


2.c 










:: 


ZnTe 


43m 


10.6 


440 


3.6 










2 


CdTe 


43m 


10.6 


800 


7 










2 


GaP 


43m 


1.06 


525 


1.3 










4 


GaP 


43m 


1.06 


255 


0.6 










3 


GaAa 


43m 


1.00 


1,500 


1 










-1 


GaAs 


43m 


10.6 


1,760 


3.7 










2 


InAs 


43m 


10.6 


2,000 


3.2 










2 


KII0PO4 


42m 


1.06 


3 


3.6 


dm 

3 


A 

3.6 







5 


KD.PO4 


42m 


1.06 


2.7 


3.2 


2.7 


3.2 






5 


KH 2 As0 4 


42m 


1.06 


3.4 


2.6 


3.2 


2.5 






4 


NH 4 H.P0 4 


42m 


1.06 


2.9 


3.15 


3. 


3.15 






5 




dm 


A 


dm 


A 


dm 


A 




ZnO 


Bram 


1.06 


43 


3.3 


13 


1.1 


14 


1.1 


4 


ZnS 


6mm 


1.06 


84 


1.9 










3 


ZnS 


6mm 


10.6 


180 


4.9 


90 


2.45 


102 


2.7 


2 


CdS 


6mm 


1.06 


186 


3.2 


96 


1.6 


105 


1.8 


2 


CdS 


6mm 


10.6 


210 


5.4 


126 


3.3 


138 


3.6 


2 


CdSe 


6mm 


1.06 


500 


3.4 










3 


CdSe 


6mm 


10.6 


260 


4.S 


136 


2.4 


148 


2.6 


2 


BaTiO, 


4mm 


1.06 


42 


1.0 


111 


2.45 


105 
dm 


2.35 

A 


5 


LiNbO, 


3m 


1.06 


250 


9 


36 


1.1 


19 


0.55 


6 


LiNbO, 


3m 


1.152 


dm 


A 


32 


1.05 


15 


0.45 


6 


Si0 2 


32 


1.06 


2.5 


1.9 










4 


AIPO4 


32 


1.06 


2.5 


2.2 










4 


Se 


32 


10.6 


380 


2.1 










2 


Te 


32 


10.6 


25,400 


4.3 










7 



946 
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Table III — Optical Rectification Coefficients 



Units of d 10"' Mil Units of A 10 -; esu 



Material 


Class 


x„ 


An 


A 




Kef. 


ZnTe 


23m 


0.694 
1.06 


3650 
1040 


14 

5 






8 

8 








dm 


A 








KH 2 PO, 
KD 2 P0 4 
NH^HoPO, 

CdS 


42m 
42m 
42m 

6mm 


0.694 
0.694 
0.694 

0.694 


132 
d»i 

700 


3.0 

A 

7 


50 
105 

d»n 

900 


3.2 
2.9 

A 
9 




8 
8 

'.) 

8 



published s.h.g. data, 7 values from optical rectification data and 50 
from electro-optic data. Definitions and conventions are discussed in 
the appendix and a separate list of references is given for the data in 
the appendix. Probable errors vary from measurement to measurement. 
It is probably safe to say that no measurement has an accuracy better 
than ±10 percent and in many cases the probable error is greater. That 
for the s.h.g. data at 10.6 n is 30 percent and except for ADP the rectifica- 
tion data is only good to a factor of 3. In addition, a few materials have 
discordant results reported by different groups and this suggests that, 
especially in the case of crystals which are difficult to grow, the data 
should be regarded rather critically. In the case of CuCl, Sterzer, 
Blattner and Miniter 20 have constructed a modulator whose behavior 
is consistent with the higher value of the electro-optic coefficient. This 
casts some doubt on the low value for CuBr reported in conjunction 
with CuCl. In the case of the linear e-o coefficient in HMT, Heilmeyer's 13 
value d 123 = 32 X 10" 9 esu is the most recent and reliable. 

The average of all the s.h.g. data is A = 3.3 X 10" 8 esu, and only 
two coefficients d 123 in HMT and d 383 in LiNb0 3 exceed 6 X 10" 6 esu 
by more than the probable error. One coefficient d 222 in LiNb0 3 is 
unambiguously less than 1 X 10" 6 esu. The sole accurate rectification 
coefficient, A 123 in NH 4 HP0 4 , is 3 X 10" 6 esu which is remarkably 
close to the values 3.15 and 3.2 X 10~ 6 esu for s.h.g. and the electro- 
optic effect. 

Whereas the s.h.g. data have a rather compact distribution about 
3 X 10 -6 esu the linear e-o data are more straggled. The mean value 
is 2.3 X 10~ 6 esu but there are a considerable number of materials 
with A < 1 X 10~ 6 esu. The difference between the averages A. hg and 
A,„ is not due to the different materials in the two lists, it persists if we 
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Table IV — Electro-optic Coefficients 



I'nits of «/ 1()-' J esu I'nits of A 10 _B esu 



Material 


Class 


X 


dm 


A 








Kef. 


HMT = N4(CH 2 )6 


43m 


n , r > 


M 


14 








10 


UMT = N4(CI1»)6 


33iii 







2.3 








11 


HMT = N4(CH,)u 


43m 




55 


21 








12 


Bi4(Ge04)j 


43m 




22 


0.8 




( m 





13 


Sodalite 


33m 




9.5 


1.8 








14 


CuCl 


33m 




28 


0.75 








10 


CuCl 


|3m 




110 


3 








15 


CuBr 


33m 




22 


0.4 


assumed 


, — 


10 


10 


ZnS 


[3m 


no 


74 


0.9 








17 


ZnSe 


43m 


.v> 


120 


0.8 








18 


ZnTe 


43mi 


(1 «( 


440 


1.5 








19 


GaP 


33m 


o.ea 


150 


0.3 








20 


GaAs 


43m 


1 (12 


215 


0.3 








21 


NaClOa 


23 


(1 fin 


2.5 


6 








22 


KiMei(S04)i 


23 




<26 


<0.1 




t = 




23 


(NHOjMniCBOi). 


23 




4.3 


0.5 




f = 




24 


(NH4) 2 Gd 2 (S04)3 


23 




;» . i 











24 


NaVO:(Cll«COO) 3 


23 




5.3 


1 3 








25 


NaaSbSVOHtO 


23 




1(1 


2 




' 




20 


1 ren Chloride 


23 




9 5 


2 7 








27 




(1mm 


(1 (13 


'/m 


A 


''11s 


A 






'/jUO 


50 


1.5 


26 


8 


)«/ll. 


28 


Zn> 


(iinin 


0.83 


07 


0.9 


34 


(1 45 


}- neK. 28 
Join -is 


(Jda 


Um id 


0.63 


110 


1 2 


48 


0.55 


Material 


('lass 


X 


urn 


A 


</.v;i 


A 


Kef. 


KUiPOi 


42in 


(1 55 


+65 


4.0 


-50 


1 7 


constant stress 


29 


KlhP(>4 


42m 




60 


3.7 








30 


KD1PO4 


42m 




160 


4.0 








31 


KHsAbOj 


42m 




77 


3 7 


81 


1 7 


constant stress 


32 


RbHtAaO* 


42m 




92 


3.5 






constant stress 


32 


NH4H2PO4 


12 in 




+55 


4.4 


-140 


3 4 


constant stress 


32.29 


NH4H1PO4 


42iu 




3(1 


3 2 






constant strain 


29 






















</iu 










</iu 


A 




A 






A 


- - 
























dm 




Ha Til ).. 


4mm 


63 


1000 


1 






320 




0.3 


+ 


33 












3.3X10* 


1.4 


constant strain 




34 












6.6X10' 


1.9 


constant stress 




35 
























dm 










./ill 


A 


(/»■ 


A 


rim 


A 


dm 


A 


'/ill 




LiNbOi 


3m 


63 


860 


4.3 


840 


2.5 


280 


1.2 


110 


0.3 


+ 


36 




1/111 


A 








SiOi 


32 


5 


3.2 


0.9 








22 


KiSjO. 


32 


) 55 


1.4 


0.4 








11 


.SrS;0 6 - HlO 


32 


) 55 


0.05 




. 






11 


CeHuOeNaHr-lljU 


32 


1 55 


65 




>assume 






11 


C8C4H4O. 


32 


0.55 


7 


1.4 


1 






11 






</jl! 


A 






C/ClTsOHh 


4 




9.7 


7 




37 
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eliminate all materials not common to both lists and may therefore, be 
either a real effect or a systematic error. 

A few materials [e.g., SrS 2 6 -H 2 0, C a H 12 9 NaBrH 2 and K 2 Mg 2 
(S0 4 ) 3 ] have very low values of A. The latter is especially interesting 
since the isomorphous (NH 4 ) 2 Cd 2 (S0 4 ) 3 and (NH 4 ) 2 Mn 2 (S0 4 ) 3 salts 
have somewhat larger values. The ammonium cadmium salt is known 
to be ferro-electric at very low temperatures and the ammonium man- 
ganese salt is also suspected of ferro-electricity (Jona and Shirane 28 ). 
More significant, perhaps, is the fact that the divalent ions have very 
nearly regular octahedral coordination (Zemann and Zemann 30 ) and 
so form a unit with near inversion symmetry and contribute little to d. 
The main contribution comes from the monovalent ions and their 
irregularly placed neighbors. The difference between the potassium and 
ammonium salts would then be due to the difference in the polariza- 
bility of the two ions. For K + the refractivity is 2.45 ccs and for NH + , 
it is 4.05 ccs (see Le Fevre 12 ). If this enters d as a cube the expected 
ratio of the d coefficients would be 4.5. The observed value is greater 
than about 5. Note also that NH 4 itself lacks inversion symmetry. 

The tabulated values of A show that Miller's rule is an excellent rough 
guide to the probable value of d. If the component is allowed by sym- 
metry 

| d-f: |:^3X 10- fl X"-X?,.A7. esu. (136) 

However, the rule by itself is not infallible. Occasionally, a value of 
d much higher than that predicted by (136) will occur. In some cases, 
(e.g., d 333 in LiNb0 3 ) this is accompanied by a very low value of another 
coefficient and it is then plausible to assume that this is due to a partic- 
ularly critical geometric configuration. In other cases, (e.g., HMT) it 
is quite clearly due to the coincidence of a number of favorable factors. 
The atoms, the molecule and the crystal all have the same symmetry 
and moreover, as we saw in Section 1, all the separate contributions to 
d have the same sign. Thus, it is likely that the value of A « 15 X 10"° 
esu for HMT represents something of an upper limit to what is possible. 

More often (136) will overestimate d. This is especially likely to 
occur if the molecules themselves, or large sections of the molecule have 
near inversion symmetry, but it may also occur if the crystal structure 
itself departs only very slightly from a centro-symmetric structure. 

VII. CONCLUSION 

If reasonably good ground state wave functions are available, the 
direct perturbation method of Section II seems most suitable as a 
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basis for calculating the magnitudes of the coefficients. It gives the 
intrinsic nonlinear coefficient 

JO '*X ' ■)*■ /10~\ 

dak = (137) 

a„e ii v 

in terms of the intrinsic low-frequency limit of the optical suscepti- 
bility x and a cubic moment in the ground state. If the electrons are 
uncorrclated, this can be replaced by 

r, a = f£ «,,. (138) 

and t iik can be obtained from the charge distribution. From (13S) we 
obtain the reduced tensor 

A, y , = ~-hh (139) 

and we can then incorporate this directly in Miller's rule using the 
observed susceptibilities x," , etc. to obtain d°f k y . 

This continuation of the basic perturbation calculation with Miller's 
rule appears to be the most straight-forward approach to the coefficients. 
Apart from the cubic moment t,, k it involves only experimental 
quantities. 

The analogy with the classical anharmonic oscillator established in 
Section IV seems most likely to be fruitful in qualitative discussions 
of the general behavior of the coefficients. It appears to have both 
empirical and theoretical justification. 

Obviously, on this basis further generalizations of Miller's rule are 
possible. For example, we might expect the fourth rank tensor d"f k \ 6 
which describes induced second harmonic generation, the Kerr effect, 
etc. to satisfy a relation of the form 

dgS 1 = A7,,Y?,-A7 t XJ,A,,.,, . (HO) 

A calculation based on fourth-order perturbation theory and a lavish 
use of sum rules leads to 

A«««1.5X lO^^fe, (141) 

e \r / 

where (r 2 ) is the mean square radius of the charge distribution and 
Qnki is the semi-invariant 

Quki = (XiX f x k x,) - 2(x i x i )(x k x l ) (142) 

if we assume that all odd moments vanish. 
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If we take (r 2 )* as 1 A this gives 

A„„«3X Hr 10 feesu. (143) 

We have seen that in the lower-order processes T/r 3 is of the order of 
2 X 10" 2 . This does not imply that Q/r 4 is of the order (2 X 10~ 2 ) 4/a ~ 
5 X 10" 3 for whereas T is nonzero only because of asymmetric molecular 
and intermodular forces, Q is nonzero even for free atoms or ions. For 
example, in the hydrogen atom Q,,-,-,/<r 2 ) 2 = 5/18 and Qua/ir 2 ) 2 = 
1/9 so that we expect A to be of the order of 3 X 10" n to 10 esu. 
For calcite with Xopuo..i = 01 and Xa. = 0.55 this gives a value of d 
between 3 X 10" 14 and 10" 13 esu. Bjorkholm and Siegman 31 have meas- 
ured 3 X 10" 14 esu. 

We have seen that the reduced tensor A,,* is proportional to the cubic 
moment, 

Talc = {XiXjX/j 

and it is therefore clearly symmetric in all its indices. This is in agree- 
ment with Kleinman's 2 hypothesis and follows from the origin of the 
nonlinear behavior in the electronic motion. 

Finally, we remark that nothing increases d like large values of the 
linear susceptibilities yet, although, the values of most allowed reduced 
tensor components A,-,-* are near 3 X 10" 6 esu they can vary by a factor 
100:1. The molecular geometry will often indicate which end of the 
range is likely to apply to a particular material. 
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APPENDIX 

We have throughout adopted a notation, originally introduced by 
Bloembergen 15 and his colleagues, in which two fields with complex 
time dependence JEfe**' and E y k e iyi produce a polarization P°e ,at at 
the algebraic sum frequency a = j8 + 7 according to 

Py° l = d^StBTfi iV * 7U . (144) 
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If the actual fields vary as cos co,/. and cos u) 2 t there will be terms in 
P at co, -f- wo and co, — w 2 obtained from (144) by letting /3 = co, , d>, , 
etc. where di, = — co, . 

This notation has several advantages in theoretical calculations for 
much the same reason that the use of complex numbers simplifies ac 
circuit theory, and for much the same reason it has a number of dis- 
advantages in calculating numerical values. For this reason, it has not 
gained general acceptance by experimentalists who tend to use a number 
of different notations, some of which, especially in electro-optics, are 
of respectable antiquity. The difference between the two notations 
introduces various factors of 2. These are independent of the general 
reluctance of physicists to state unequivocally whether they are using 
peak or rms fields. In particularly fertile ground, these various factors 
can luxuriate and blossom as factors of 8 in the final answer. We use 
peak fields in all our definitions. 

If the applied field is 

F(l) = (0, F, cosco/, /'';, cos Qt), (145) 

it has components E\ = E\ = hF 2 , El = E" = %F a and the 1 com- 
ponent of P.(t) is 

P,(0 = il^W"" +<*} 4- |{<OVV' Cu+D) ' + cc\ 

+ ttd:JF 2 F :t e i( - QU + ec) + lMlF 3 F s e Ua - au + cc\ , (146) 

where we have used d'{tl = d£j an d suppressed the first superscript, 
which is always the algebraic sum of the second and third superscripts. 
We can also write (146) as 

p.(0 = IWS + rfSlWi cos (« + «)' 

+ }{#» + d,f 2 \F 2 F, cos (co - G)/. (147) 
If co = 12, this gives 

P.(0 = h\d\tr + d\°r\F 2 F, cos2to/ + *{<*?£■ + d\%:)F 2 F* . (148) 
Now the usual experimental definition would be 

Pi(t) = (<*& + d*;. 2 )F 2 F, cos 2co/ + (d° l2 , + d° ia2 )F s F a (149) 
and so we see that 

s.h.g. dj- = \d\ir (150) 

rectification d" llk = &ft£ . (151) 
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If we let ft = 0, we have 

P.(/) = (« - + dST)JW coso,/. (152) 

The experimental definition reads 

Pl(t) = <&|7|F| COS CO/ (153) 

i.e., 

5 Xl2 = d" 23 F, , (154) 

so that it is possible to contract the last two suffices according to the 
scheme 

U_»l 22 ->2 33->3, 32 = 23->4, 31 = 13->5, 12 = 21-»6. 

(155) 
Thus, d% (i = 1 • • 3, p = 1 • • • 6) represents djj[ and, for example, 

"25 — "231 = "213 • 

P, cos 2cu/ = (d?2 3 + d 2 32 )F 2 F a cos 2co/ 

= 2d 2 i23 F 2 F 3 cos 2co/ = cti a t -2F 2 F 3 cos 2u>/. (156) 
It is therefore common to define the "vector" 5 

5 - ff, , ff 2 , ff, , ff 4 , S 3 , CF 6 = F? , F 2 2 , Fl , 2«, , 2 V, , 2F,F 2 (157) 
so that 

P f = E <W, . (158) 

i>=i 

With this notation d% = d?£ = d$& ^ d? 2 % + dun Dut als0 ^ = d i" • 
In the electro-optic case, the subscripts referring to optical fields 
can be contracted 

aiT,; - = dt ik = C p = 1 ■ ■ ■ 6, fc - 1.3. (159) 

The alternative ordering (155) leads to d£, . 

Note that in this case, since d operates on two distinct fields, one 
optical the other dc, there is no possibility of constructing a "vector" 
such as 3\ The sum implied in the definition of d° k is 

5x P = &Xu = td;X c . (160) 

Electro-optic data are often presented as coefficients r pk in the sus- 
ceptibility ellipsoid. If n is the refractive index (assumed isotropic), 
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4?r ,„. 4t _,„ 8tt Illl(J o ,.„«•, 

'/-*• = -"~4«m- = ——iO-iik = ~—\d iik . (161) 

The dimensions of d and r are those of an inverse field. 

In the MKS system, the units are meters per volt, in the cgs system 
they are centimeters per stat-volt. One MKS unit is 3 X 10 4 esu and 
so numerical values of d in esu are the larger numbers. 

We have not discussed the influence of a mixed use of rms and peak 
fields but we note that if rms fields are used throughout the values of 
the coefficients will all be V2 times larger than if peak fields are used 
throughout. No one is, however, likely to use an rms dc field. 

Experimental values of the electro-optic coefficients are usually 
expressed in absolute units and the only ambiguity that can occur is 
associated with whether the measurements were made at constant stress 
(undamped) or constant strain (clamped). It is safe to assume that 
constant stress is implied by the absence of any definite statement to 
the contrary. 

Second harmonic coefficients are sometimes given in absolute units 
but more often relative to the coefficient d^ in KH 2 P0 4 . An absolute 
measurement of this by Ashkin, Boyd, and Diedzic 32 gave 

dloy = \d\tr = 3 X 10 ~° esu, 

but this is now believed to be too large. The most recent measurements, 
Francois 5 , Bjorkholm,™ give 

dltx = hdl™* = 1.38 X 10 " 9 esu =fc 12 percent 

for the coefficient in NH 4 H 2 P0 4 . Relative measurements show that it is 
identical in KDP and ADP. We have used a rounded off, compromise 
value 

KDP4", = yin* = 1.5 X 10"" esu (162) 

in compiling the tables. It affects all values of d lw at optical frequencies 
but not at 10.6 p. 

It will be apparent that in comparing theory or experiment with 
experiment, considerable care is needed to be sure that like definitions 
are being compared with like. 
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